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Abstract. A class of non-canonical inflationary models is identified, where the 
leading-order contribution to the non-Gaussianity of the curvature perturbation is 
determined by the sound speed of the fluctuations in the inflaton field. Included in this 
class of models is the effective action for multiple coincident branes in the finite n limit. 
The action for this configuration is determined using a powerful iterative technique, 
based upon the fundamental representation of SU(2). In principle the upper bounds 
on the tensor-scalar ratio that arise in the standard, single-brane DBI inflationary 
scenario can be relaxed in such multi-brane configurations if a large and detectable 
non-Gaussianity is generated. Moreover models with a small number of coincident 
branes could generate a gravitational wave background that will be observable to future 
experiments. 
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1. Introduction 

The quest to realise inflation within string/M-theory continues to attract considerable 
attention. The Dirac-Born-Infeld (DBI) scenario of the compactified type IIB theory is a 
well-motivated model, in which inflation is driven by one or more D-branes propagating 
in a warped 'throat' background [HEIEIIUEIEICOIEIEI ED] • Such a background is 
generated by the non-trivial form-field fluxes over the internal dimensions. (For recent 
reviews, see [HI [12j E3 HU EE1 EEl EE]-) In the simplest version of the scenario, the 
inflaton parametrizes the radial position in the throat of a single D3-brane. The brane 
dynamics are determined by the DBI action in such a way that the inflaton's kinetic 
energy is bounded from above by the warped brane tension. The regime where this 
bound is nearly saturated is known as the 'relativistic' limit. 

Recently relativistic DBI inflation has come under considerable pressure when 
confronted with cosmological observations. Baumann & McAllister (BM) and Lidsey 
& Huston (LH) have shown that the ratio of the amplitudes of the tensor and scalar 
perturbations generated during inflation is bounded from above by r<lCT' [18, 19J. 
However in ultra-violet (UV) versions of the scenario, where the D-brane is moving 
towards the tip of the throat, the tensor-scalar ratio is also bounded from below, 
r>0.1(l — n s ), where n s denotes the spectral index of the scalar perturbation spectrum 
[19J. The two bounds on r are incompatible if n s ~ 0.95, as currently favoured by 
Cosmic Microwave Background (CMB) observations [20] |2T] . 

The purpose of the present paper is to investigate whether the upper bounds on r 
can be relaxed in more general DBI inflationary scenarios. A natural extension to the 
single D3-brane model is to consider a Dp-brane wrapped around a (p — 3)-cycle of the 
internal space. For example, Becker et al |22j have proposed a model where inflation is 
driven by a D5-brane. In this case, the range of allowed values for the inflaton becomes 
independent of the throat charge, N, which weakens the upper bound on the tensor- 
scalar ratio to r<0.04. Strictly speaking this is only true for the D7-brane case, since 
the wrapped D5-brane imposes A0 ~ iV" 1 / 4 . However in arriving at this bound, it 
was assumed that the backreaction effects of any fluxes in the throat were negligible. 
Kobayashi et al [23] considered both D5- and D7-brane models, but concluded that the 
former case required an excessively large background charge in order to relax the bounds 
on r. Whilst this is highly constraining, it is still much better than the case for single 
D3-branes which require a much higher background charge - and therefore are effectively 
ruled out as a predictive model. Thus wrapped brane configurations are preferable to 
single brane models. However the difficulty in these models is that the backreaction is 
no longer under control. 

Alternative ways to relax these bounds have been proposed, including theories based 
upon multi-field models [21], the addition of angular momentum as another degree of 
freedom [25] and using different throat geometries [26]. However it must be noted that 
the extra degrees of freedom introduced in these models do not solve the problem. The 
bounds are relaxed only by a small fraction, and therefore these models should still be 
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regarded as being unsatisfactory since they require an extreme amount of fine tuning in 
order to work. 

Another alternative possibility is to consider multiple brane configurations^!- In the 
case where n branes are localised initially at equal distances I > l s and subsequently 
follow the same trajectory, the effective theory is equivalent to that of n copies of the 
action for a single brane. A more general initial condition, particularly for branes 
created in the infra-red (IR) region of the throat [H [28l [29] , is that the branes should 
be separated over a range of scales, with a subset being coincident and the remainder 
being widely separated. 

Our approach in this paper is two-fold. We begin in Sections [2] and [3] by noting 
that the upper bounds on the tensor-scalar ratio arise due to the special algebraic 
properties of the DBI action. We then adopt a phenomenological approach in Section [4] 
and identify a general class of non-canonical inflationary models, where the leading-order 
contribution to the non-Gaussianity of the curvature perturbation is determined entirely 
by the speed of sound of the inflaton fluctuations. In these models, the bounds on r can 
be relaxed if significant non-Gaussianities are generated. This class of models includes 
the relativistic limit of the action for n coincident D3-branes, which originates from a 
UV complete theory. This motivates us to develop the theory of multiple coincident 
branes further in Section [51 We find that for a finite number of branes, the effective 
action for n coincident branes can be derived and the backreaction kept firmly under 
control. In Section [6] we find that such models can in principle lead to a detectable 
gravitational wave background if the number of coincident branes is sufficiently small. 

Units are chosen such that H = c = 1 and M P = (8nG)- 1/2 = 2.4 x 1(T 18 GeV 
denotes the reduced Planck mass. 



2. Non-Canonical Inflation 



The low-energy, world-volume dynamics of a D3-brane in a warped background is 
determined by an effective action of the form 



S 



d xi 



\9\ 



M 2 



(1) 



where R is the Ricci curvature scalar, X = — \g^ y V ^(pV u 4> denotes the kinetic energy 
of the inflaton field 0, and the function P((p,X) is referred to as the 'kinetic function'. 

We assume that the four-dimensional universe is spatially flat and isotropic and 
sourced by an homogeneous inflaton field, <p = <p(t), with energy density E = 2XP t x—P, 
where a subscripted comma denotes partial differentiation. We further assume that 
the inflaton dynamics generates a quasi-exponential expansion of the universe, where 
e = -H/H 2 < 1. 



% In certain limits this approach is actually dual to considering wrapped branes [27 . 
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It proves convenient to define two parameters in terms of the kinetic function P 
and its derivatives [301 



r 2 = ( 2 ) 

° s ~ P x + 2XP XX ' [Z) 

The first parameter, c s , determines the sound speed of fluctuations in the inflaton field. 
This can be significantly less than unity, in contrast to slow-roll inflation driven by a 
canonical field such that P yX = 1- 

The amplitudes of the scalar and tensor perturbations generated during inflation 
are given by [32] 

respectively, and the ratio of these amplitudes is defined as (32] 
Vl 

r = ^ = lQc s e. (6) 
' s 

The WMAP3 normalization of the CMB power spectrum implies that V\ = 2.5 x 10~ 9 
and the experimental upper bound on the tensor-scalar ratio is r < 0.55 [20] . 

Deviations from Gaussian statistics in the curvature perturbation, TZ, are 
parametrized in terms of the non-linearity parameter, /nl, which is defined by TZ = 
TIg + §/nl(7^g — (^g))j where the quadratic component represents a convolution and 
TZg denotes the Gaussian contribution [33]. In the limit where the three momenta 
have equal magnitude (corresponding to the equilateral triangle limit), the leading-order 
contribution to the non-linearity parameter is given by [5H EI] 

One should note that the sign convention is that employed by the WMAP data set. Data 
from WMAP3 imposes the bound |/nl| < 300 on this parameter [20]. The corresponding 
bounds for other triangle configurations may be much tighter than this and this may be 
particularly relevant if non-Gaussian signatures have indeed been detected in the CMB 
[351 [36] • The more recent WMAP5 data set [21] improves on this bound somewhat, and 
also indicates that it is distinctly asymmetric. At the 95% confidence level, the bound 
on the equilateral triangle becomes —151 < /nl < 253. 

Eqs. (jlj) and (jSJ) imply that the variation of the inflaton field during inflation is 
related to the tensor-scalar amplitude by [371 [IB] 
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where M = J dt H denotes the number of e-foldings. We will refer to the epoch of 
inflation that can be directly constrained by cosmological observations as 'observable 
inflation' and will assume that this phase occurred when the brane was located within 
a throat regional Observable inflation corresponds to no more than about 4 e-foldings 
of inflationary expansion, AjV* — 4. The total variation in the inflaton field between 
the epoch of observable inflation and the end of inflation is then given by 



A0 ; - / - ^ 



M P V 8c,P 



.X 



Ks , (9) 



r P V \V2 rKnd / „ \ 1/2 
'"xi X 



where 

X feJ ^ (10) 

If r/(c s P t x) varies sufficiently slowly during observable inflation, the corresponding 
change in the value of the inflaton field is given approximately by [371 [T8] 

A0\ 2 (AAQ 2 



Mr J S • ./'v.' ' l L1> 

3. Theoretic Upper Bounds on the Tensor-Scalar Ratio 

The ten-dimensional metric of the warped deformed conifold inside a throat region has 
the form 

ds\ = h 2 (p)dsl + h- 2 (p) (dp 2 + p 2 ds 2 X5 ) , (12) 

where the 'warp factor' h(p) is a function of the radial coordinate p along the throat 
and X 5 denotes a five- dimensional, Sasaki-Einstein manifold. In many scenarios, the 
ten-dimensional manifold (fl2"|) can be approximated by the product Ad £5 x X 5 , where 
AdSc, represents five-dimensional, anti-de Sitter space. In this case, the warp factor is 
given by 

L Vol(X 5 )mf 1 ; 

where L denotes the AdS^ radius of curvature, Vol(X 5 ) is the volume of the five-manifold 
X5 with unit radius, iV is the D3-brane charge in the throat, g s is the string coupling 
and m s is the string mass scale. Generally the value of the inflaton field is determined 
by the radial position of the D3-brane in the throat, (p = py/Ts, where T 3 = m 4 /[(27r) 3 ^ s ] 
is the brane tension. 

The four-dimensional Planck mass is related to the volume of the compactified 
Calabi-Yau three-fold, V e , such that Mp = Vqk^q, where nf = \(2ti) 7 g 2 s / ml = n/T 2 for 
a D3-brane [JJ. Hence the Planck mass is bounded from below by the volume of a throat 
region, Mp > V^th^io 2 ' where V^th^^e denotes the throat volume. For an AdS 5 x X< 



§ We denote the values of all parameters evaluated during observable inflation by a subscript '*'. 
I We parameterise the Planck scale in terms of the D3-brane tension out of convenience, and note that 
there is no physical relationship between the two. 
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throat, Baumann & McAllister (BM) exploited this inequality to derive an upper limit 
on the maximum variation of the inflaton field in the throat, A0 max < 2Mp/^/N, which 
leads to the corresponding limit |A0|* < 2Mp/\HV [15] . Combining this with the 
constraint (ITT]) therefore yields an upper limit on the tensor-scalar ratio: 
32 



n< mf s {CsPx) * ■ (14) 

Two of the authors (Lidsey & Huston, LH) derived a complementary bound on the 
tensor-scalar ratio by noting that during observable inflation the brane spans a fraction 
of the throat volume [19] 



|AV 8 ,,| ~ Vol(X 5 )l^M (15) 
and, since |AVe,*| < ^6,th, it follows that 

a<a 2 r 3 ^ (A P ,) 2 

It was then assumed that the fractional change in the value of the inflaton field during 
observable inflation was less than unity: 

|A0*|<0*. (17) 

This condition is necessarily satisfied in UV versions of the scenario, where the brane is 
moving towards the tip of the throat, but must be assumed as a further constraint in 
IR versions where the brane moves out of the throat, since in these latter cases 0* could 
be very small. Combining the limit (fI71) with Eq. ([15]) then implies that 

|AV 6 ,,|>Vol(X 5 )^^ (18) 

and substituting this constraint into the bound ([16]) yields the condition 

^v <T7 izw^y. (19) 



r.<^rAJ^X M- (20) 



MpJ, Vo\(X 5 )\M P y 
Finally substituting the constraint (ITT]) yields the upper limit [TJJ 

1/3 4/3 

(AAO 2 VVol(X 5 ); \M P 
Comparison of the limits ( ITTj) and ( 120]) implies that the LH bound is the stronger 
of the two when 

'^< 2 °™*) 1/3 (^)" /3 ^- < 21 > 

For typical field-theoretic values Vol(X 5 ) ~ 0(n 3 ), m s ~ 0.1M P and g s ~ 10~ 2 , this 
implies 

hf'N < 300^i . (22) 

- /V eff 

In the following Section, we identify a class of models in which these bounds could 
be relaxed. 
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4. Relaxing the Upper Bounds on the Tensor-Scalar Ratio 

In the standard DBI scenario, the kinetic function defined in Eq. (OQ) takes the form 



P(0,X) = -T(0)v/l-2T-i(0)X + T(0) - V {<j>) , (23) 

where T(</>) = T 3 /i 4 (0) is the warped brane tension and V(0) is the inflaton potential. 
Typically in warped compactifications of IIB supergravity, this potential is determined 
by the relevant fluxes and brane interaction terms. We will ignore the precise origin and 
form of this potential, but simply note that it is highly sensitive to the string theoretic 
construction. For the purpose of this note we will simply treat it as an arbitrary function 
of the inflaton field. (See, for example, Ref. [9] for a discussion on the precise form that 
the inflaton potential may take.) 

The standard DBI scenario (|23|) is algebraically special, in the sense that the kinetic 
function satisfies the constraints 

c s P, x = 1, A = ~ (\ - 1 ) . (24) 



2 W , 

It follows that the bounds (|14p and (T2"Uj) on the tensor-scalar ratio could in principle be 
significantly relaxed in models where (c s Px)* 3> 1- In view of the second relation in 
Eq. (1241) . it is of interest to take a phenomenological approach and consider the more 
general class of models where 

\ - 1 = aA , (25) 

for some positive constant a. Moreover since a large non-Gaussian signature in the 
curvature perturbation is typically generated in models where the sound speed of 
fluctuations is small, we will begin by considering scenarios where the kinetic function 
satisfies the inequalities: 

X 2 P tXXX > XP )XX > P x . (26) 

In these limits the constraint (j25l) reduces to the third-order, non-linear, partial 
differential equation 

P xx = t:P,xP,xxx ■ (27) 
o 

Changing the dependent variable to Q = P, xx / P }X reduces Eq. ( 1271) to 

aQ, x = {6-a)Q 2 , (28) 

and it is straightforward to integrate Eq. (|28]1 exactly. The remaining integrations can 
also be performed analytically and the general solution to Eq. (1271) for a ^ 6 is given 

PfaX) = fM [1 - f2(<P)X] m - / 3 (0) , (29) 

% The special case a = 6 results in an exponential dependence of the kinetic function on X, and is 
therefore an example of a higher-derivative theory. However we do not consider this model further, 
since it does not lead to a weakening of the gravitational wave constraints. 
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where fi(4>) are arbitrary functions of the scalar field and 

2(a-3) , , 

m = ■ 30 

a — 6 

It can be verified that the inequalities fl26|) are satisfied in the 'relativistic' limit, where 
X ~ 1//2- We consider the inflationary dynamics in this limit in what follows. For 
completeness, we note that Eq. (125|) can be solved analytically in full generality and 
the solution is presented in the Appendix. 

The standard DBI scenario is recovered for m = 1/2. More generally however, Eq. 
f l29|) implies that 

m flf2 n n v \(2m-l)/2 

c s Px ^ , (1 - f2Xy , (31) 

a/2(1 - m) 

2(1 — m) 

when X ~ I//2. Self-consistency therefore requires m < 1. Moreover we find from Eq. 
(ED that 

f a 5 (59 ~ 55m) 

/NL "l— frX> P= 486 ' (33) 
a 5 / 59 — 55m \ , , 

~ ~cf' ^ 972 (^rJ- 

Hence substituting Eqs. (EU) and ([33]) into the BM bound flU]) and the LH bound (J2Uj) 
implies that 

r^J^^ 112 (-^Y' 2m)/2 (35) 



^ / T 3 \ 1/3 4/3 {-m)hf 2 (1 ' 2m)/2 

r * < (AAf)l \Vo\(X 5 )J \M P J Ml-m) V P J { ' 



and 



respectively. 

We conclude, therefore, that the upper limit on the tensor-scalar ratio could be 
significantly relaxed if m < 1/2, since the non-linearity parameter is at present only 
weakly constrained at /nl > —151. Although it is possible to phenomenologically 
construct a model which has a value of m in this range, it is clearly preferable to 
identify UV complete models that satisfy this requirement within a string theory context. 
Unfortunately this is quite difficult to achieve since the inflaton will either be associated 
with an open or closed string mode. The open strings are governed by relativistic 
actions of the DBI form, whilst closed strings arise from compactification of Einstein 
gravity and are typically put into canonical form. However there do exist classes of 
open string models which satisfy the above requirement, namely those associated with 
multiple coincident branes. 

More specifically, if the branes are spatially separated, the effective action is 
algebraically equivalent to that of a single brane. It will therefore not satisfy the bound 
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on mEI- Similarly it was shown in [38] that n coincident branes, in the large n limit, will 
also fall into this class of models. On the other hand, if it is assumed that n is finite, 
the special properties associated with the matrix degrees of freedom become important 
and this results in a kinetic function satisfying m < 1/2. We will discuss this in more 
detail in the following Section. 



5. Action for Multiple Coincident Branes 

We have seen how the form of the kinetic function P can significantly change the strength 
of the LH bound on the tensor-scalar ratio, depending on its explicit form. One model 
in which a suitable form for P is realised is the multiple coincident brane model as 
outlined by Thomas & Ward [38] . 

The world-volume theory for coincident branes is not fully known, although 
a number of proposals have been made. We will restrict our analysis to Myers' 
prescription, since this has been extensively discussed in the literaturqj [iQj E]. In 
general the open string degrees of freedom for n coincident branes combine to fill 
out representations of U{n) (as opposed to U(l) n in the case of separated branes). 
This introduces a non-Abelian structure into the theory. In the single brane case, the 
fluctuations of the brane are characterised by induced scalar fields on the world-volume. 
However for multiple branes these scalars must be promoted to matrix representations 
of some gauge group. 

Typically the transverse space of any given compactification will always admit an 
SO (3) isometry. We can therefore choose our scalars to transform under representations 
of the algebra of SO (3) ~ SU(2) by making the identifications 

0* = Ra i i = 1, . . . 3 , (37) 

where R is some scale with canonical mass dimension, and the a % are specified to be the 
irreducible generators satisfying the commutator 

2ie ijk a k , (38) 



\a , a 



and the conditions 

-Tr(aW) = CS ij = {n 2 - l)8 ij , (39) 
n 

where C is the quadratic Casimir of the gauge group. The irreducibility condition 
corresponds to the configuration being in the lowest energy state. It is therefore an 
additional fine-tuning of the initial conditions. 

The Myers prescription requires a symmetrised trace (denoted STr) to be made 
over the gauge group. This implies that the symmetric averaging must be taken over 
all the group dependence before taking the trace. For n>l, the symmetric trace can 
be approximated with a trace, which results in the usual DBI action multiplied by a 

+ In this discussion, we are ignoring the non-trivial backreaction of these branes on the background, 

and therefore one should be careful about the range of validity of the effective action. 

* There is also a proposal by Tseytlin [39] for the non-Abelian theory of coincident D-branes. 
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potential term (as described in [3SJ[2n]). However for finite n, the symmetrisation clearly 
becomes more important and it is essential that we have some means of performing this 
operation. Recently a prescription for the symmetric trace at finite n was proposed 
[121113], using highest weight methods and chord diagrams. 

The result is that the STr acts on different spin representations of SU(2) in the 
following manner: 

n/2 

STr^oty = 2(2g+ l)J^(2i- l) 2 ", n even , (40) 

i=i 

(n-l)/2 

STriaWy = 2(2q + 1) ( 2i f^ nodd - ( 41 ) 

i=i 

In order for the solution to converge in this prescription, it is also necessary to 
modify the definition of the radius of the SU(2) sphere. In the large n limit, this is 
given by 

p 2 = \ 2 R 2 -Tr(a i a i ) = \ 2 R 2 C , (42) 
n 

where A = 27rZ 2 = 2nm~ 2 , whereas for finite n, it becomes 

" 2 =" 2r2l ^A^S^)^ x2r> - i)2 - (43) 

This converges to the large n result in the appropriate limit. This point is important, 
since the warp factor of the four- dimensional theory is typically of the form h = h(p). 

The resulting kinetic function for n coincident branes in the finite n limit is therefore 
given by 

P = —T 3 STr ^ 4 (p) (-ZR^Y^aW)^^ 2 ^ + V(p) - h\p) J , (44) 

where 

Z S Cf) ^ r(3/2 ™ (1 + g) - («) 

Note that the second and third terms in Eq. ( ]44"1) are singlets under the STr and 
therefore contribute terms proportional to n. The physics of these branes away from 
the large n limit is particularly interesting as discussed further in [381 127] . 

The simplest case to consider is that of two coincident branes. However the form of 
the STr prescription implies that all other solutions for n > 2 can be deduced entirely 
from the n = 2 solution by a recursion relation. In order to see this, let us define 

1 + 2Y - (2 + 3Y)ZR 2 ^ 



2T,h 




VTTyVi - ZR 2 

2Y-YZR 2 ) \ 



vTT7(l - ZR 2 ) 3 / 2 1 



(46) 
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These quantities correspond to the pressure and energy density functions when n = 2 
which arise solely from the DBI sector of the action. The full pressure and energy 
densities are then given by P = P2 — 2T 3 {V — h 4: ) and E = E2 + 2T 3 {V — /i 4 ), respectively. 
Since the symmetrised trace acts differently on the differing spin representations of 
SU(2), we should expect this structure to follow through in the recursion relation. 
Indeed, we find that for odd n 

/(n-l)/2 \ 

P n° )= I E P2mfZ,{2kfY] I -nT 3 (V-h*), 

/(n-l)/2 \ 

E n )= ( E E 2 [(2k) 2 Z,(2k) 2 Y] I +nT 3 (V-h*), (47) 
whilst for even n we find that 

/ n/2 \ 

P n E) = l£M2*-l) 2 £,(2*-l) 2 Y] J -nT 3 (F-/> 4 ), 

/ n/2 \ 

^ = ( E - ^ ( 2fc - ^ I + - h ') ■ ( 48 ) 

For example, we can employ these recursion relations to obtain the solutions for 
n = 3: 

p= _ 2T J ni+g-gzgft+gTA _ _ 

V VI - 4Zi?Vl + 4F / 

which agrees precisely with the result computed by direct expansion of the STr 
prescription. Furthermore the n = 4 case is given by 

, h 4 (l + 2Y -ZR 2 (2 + 3Y)) 

-21 a 



Vi+Yvi - ZR 2 

^(l + l^-W + ^A - 4T 3 (K - , 
VI + 9FV1 - 9ZR 2 J 

E -[ ft '" +2i '-™'' 1 + V m(v ht) (50) 

1 VTTY(i - ZR 2 f/ 2 y/T+W(i - gxi? 2 ) 3 / 2 J 

It it clear that the relevant functions increase in complexity as n increases, since there 
are progressively more terms to include in the STr expansion. However Eqs. fj47]) and 
( f48|) represent the most general solutions. 

One should also be aware that the backreaction of multiple branes will typically 
introduce corrections of the form n/N, therefore it is important for this ratio to be 
small in order for us to trust the supergravity analysis. Typically we can argue that 
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wrapped branes are dual to multi-brane configurations when we are in the limit that 
n » 1. However since we also wish to keep N » 1 we must tune the solution so 
that n/N << 1 is satisfied. Therefore the origin of the backreaction effects is much 
clearer from this perspective. One can compute the 1/N corrections to the multi-brane 
action in the large n limit [27] which, in the dual picture, correspond to backreactive 
corrections to the wrapped brane models. It would certainly be more useful to develop 
both these models in more detail. 



6. Bounds on the Tensor-Scalar Ratio for Multi-Brane Inflation 

The last terms appearing in the summations of Eqs. (jUJ and ( 1481) correspond to the 
k = (n — l)/2 term when n is odd and to the k = n/2 term when n is even. This implies 
that for all n, these terms can be expressed in the form 



h 1 



P = -2T, 



1 + 2{n - l) 2 Y - [2 + 3(n - l) 2 Y}(n - 1) 2 ZR 2 



1 + (n - 1) 2 Y\ 1 -in- 1) 2 ZR 2 



.-nT 3 (V-h A ) .(51) 



Inspection of Eqs. fj46l)-( )48l) implies that the relativistic limit is realised for any 
finite number of branes when (n — 1) 2 ZR 2 — ► 1. In this case, the dominant contribution 
to the summations appearing in Eqs. (1471) and (HBl) will arise from the last term, Eq. 
( l5Tj) . In this limit, therefore, the kinetic function appearing in the effective action 
simplifies to 

P = 2T 3 jfcVl + (n - l) 2 Y (l - ' ^-nT 3 (V-h 4 ) , (52) 

where 



Y - {n _ ijipuj [ h ) • < 53 > 

Zfi2s (n-lW Y ' (M) 
and we have effectively imposed the relativistic condition 

X ~ l -T^ , (55) 

in the numerator of Eq. ( l5Tl) . For the n = 2 and n = 3 cases, we have verified 
by direct calculation that when one calculates the speed of sound and the non- 
linearity parameter ([7J from the general expressions (14*6!) and (1471) and then imposes the 
relativistic limit (I55p . one arrives at the identical result by starting explicitly with Eq. 



At this point we should consider the validity of the function in Eq. ( 1521) . Using 
the recursion relations defined in the previous section, we see that in the large n limit 
the kinetic function converges to the corresponding function defined in the large n limit 
in [38]. This is not the same function as that for n separated branes, as the matrix 
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degrees of freedom lead to an additional potential term for the scalars. However it does 
belong to the same class of models with m = 1/2. We have verified this convergence 
numerically since the algebraic sums are unfortunately not tractable. The key point is 
that there must exist some value of n, beyond which the function appears to look more 
like the standard DBI action, rather than the approximate form proposed in (1521) . For 
a range of background solutions, the numerics suggest that the approximation is valid 
up to terms of O(10). Since there are a large number of parameters in the theory, it is 
possible to find solutions where n ^> 10. However we will then be forced to generate a 
larger background flux, which will result in a situation where even the conformal Calabi- 
Yau condition is no longer valid. In view of this, we focus on the sector of the theory 
where n < 10, which implies that the backreaction is under control and that the kinetic 
function is still of the required form. 

Eq. ( 1521) is precisely of the form given by the general solution ( 1291) . where m = — 1/2 

and 

2 



= Zr 3 h*y/l + (n-l)*Y, f 2 {<j>) = . (56) 

We may therefore immediately conclude from Eq. (|34l) that /nl — — 0.3/c 2 . Moreover, 
since f3 ~ 0.9 in this scenario, Eqs. ( 13TT) and ( 1331) reduce to 



c s P x ~ -1.3V1 + (n- 1) 2 17nl • (57) 

We first consider the LH bound (1201) . This applies at least for all UV scenarios. 
It follows after substitution of the relativistic limit (|B3|) into the scalar perturbation 
amplitude, Eq. (TJJ, that 

S 50 T 3 h^l + (n - l) 2 Y /nl ' 



Substituting the tensor-scalar ratio (JoJ) into Eq. (jooj) then results in a constraint on the 
magnitude of the warp factor during observable inflation: 

hi -1 r 2 Vl , , 

M| 2T 3V /l + (n-l) 2 F /nl ' V ' 

Eqs. ( IBTj) and ( 1391 may now be substituted into the LH bound ( 1201) to yield 

r * < (AAT)6 Vol(X 5 ) 7 5/nl • (bUj 

It is clear that the parameter Y must be sufficiently large if the tensor perturbations 
are to be non- negligible. For the AdS^ x X 5 throat, this parameter takes the constant 
value 

A-n 2 g N 

YAdS = 5T= i) 4 Voi(x 5 ) • (61) 

In what follows, we chose natural field-theoretic values for the volume, Vol(Xs) ~ 7r 3 , 
and the string coupling, g s ~ 10~ 2 , and further assume that (n — l) 2 Y ^> 1. It is possible 
that observations will probe a range of scales AM. ~ 1, but it is more realistic to require 
that the tensor-scalar ratio should not change significantly over the entire range of scales 
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that are accessible to cosmological observation, which corresponds to AA/"* ~ 4. After 
substitution of the above values, therefore, the bound fl60|) simplifies to 

r, <2.8xl0- 13 — ^-r -f* L . (62) 



Global tadpole cancellation constrains the magnitude of the background charge TV 
in terms of the topology of a Calabi-Yau four- fold such that iV < x/24, where x is the 
Euler characteristic of the four-fold jSJ HHJ 06], H7J HH1 09]. The maximal known value 
of the Euler number for such four-folds arises from hypersurfaces in weighted projective 
spaces and is given by x = 1, 820, 448 [49]. This implies the upper limit of 

N < 75852 (63) 

for known solutions, although in principle higher values are possible. Imposing the 
WMAP5 bound /nl > —151 in ( |62l) and noting that n > 2 for consistency then implies 
an absolute upper limit on the tensor-scalar ratio: 

r* <5xl(T 4 . (64) 

This limit is below the sensitivity of the Planck satellite (r>0.02) [50J. On the 
other hand, the projected sensitivity of future CMB polarization experiments indicates 
that a background of primordial gravitational waves with r*>10 -4 should be observable 
[511 [52] . In view of this, it is interesting to consider whether a detectable gravitational 
wave background could in principle be generated in this class of multi-brane inflationary 
models. We find from ( 162]) that this would require 

n < 1 - 5.3 x lO^v/iV/NL < 1 - 0.014/ NL , (65) 

where the theoretic limit (163]) for known compactifications has been imposed in the 
second inequality. We may deduce, therefore, that since we require n > 2 for consistency, 
a detectable tensor signal will require /nl < —70, which implies that an observation 
of the tensors should also be accompanied by a sufficiently large - and detectable - 
non-Gaussianity. In other words, this class of models could be ruled out if tensors 
are observed in the absence of any non-Gaussianity. On the other hand, the current 
limit of /nl > —151 implies that n < 3 is required for the tensors to be observable. 
Consequently, if tensor perturbations are detected, this would rule out all models 
with n > 4 or, alternatively, would require presently unknown configurations with N 
exceeding bound (fS3~I) . 

In the above analysis we assumed that the string coupling took the value g s ~ 10~ 2 . 
For the AdS^ xl 5 throat, the bound ( 1601 depends proportionally on g s and can therefore 
be weakened by allowing for larger values of the string coupling. For example, increasing 
this parameter by a factor of 4 to g s ~ 0.04 (so that it is still in the perturbative regime) 
relaxes the limit on the number of branes for the tensors to be detectable to n < 5. 
Similarly, considering a smaller value for the volume of the Einstein manifold X 5 will 
also weaken the upper limit. 

Let us re-iterate that this limit on n is well within the regime of validity for the 
theory, which we have argued is self-consistent for n < 10. Moreover since the constraint 
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(I65p arises using the absolute maximal bound on the known Euler characteristics, it 
suggests that in realistic scenarios n will always be much smaller than this. Indeed, one 
could argue that only the n = 2 and n = 3 theories are likely to be valid over a large 
distribution of the flux landscape. 

We must also ensure that our approximation (n — 1) 2 Y ^> 1 is valid for consistency. 
For the parameter values we have chosen this requires that g s N ^> (n — l) 2 and 
this is satisfied if the condition (165]) holds. Note also that we require N ^> n for 
the supergravity approximation to be under control and for backreaction effects to be 
negligible. This is also satisfied when (165|) holds. 

For completeness we should also consider the BM bound ( fl4l) for this class of models. 
This is given by 

42 . 

r* < - j^t V / l+F= 1 Wnl (66) 



and, in the case of an AdS§ x X 5 throat, simplifies to 

r < ^ (67) 

Comparing the limits (R)2"j) and (16T|) implies that the LH bound is stronger than the 
corresponding BM bound if 

n > 1 - 5.5 x 10- 14 AT 3 /Xff/NL (68) 

and this condition is always satisfied if 

- 5.5 x 10- 14 iV 3 / 2 jV e ff/NL < 1 • (69) 

Moreover, the bound ( 1691) will itself be satisfied for all values of /nl and N if it is 
satisfied when the limits /nl = —151 and iV = 75852 are imposed. Hence, we conclude 
that the LH bound is stronger for Af c s < 75. In general, it is difficult to quantify 
the magnitude of N e g without imposing further restrictions on the parameters of the 
models and, in particular, on the functional form of the inflaton potential. However, if 
the ratio e/P,x remains approximately constant during the final stages of inflation, one 
would anticipate that A/" e ff<60. Nevertheless, if N 75852, the bound (!68j) will only 
be violated for n < 3 if J\f e g 3> 60. 

Finally, it should be emphasized that the derivation of the LH bound underestimates 
the Planck mass by assuming that the volume of the throat is much smaller than the 
volume of the compactified Calabi-Yau three-fold. It is likely, therefore, that the actual 
constraint on r would be much stronger. Consequently, although the bound (|65|) does 
marginally allow for detectable tensors if n is sufficiently small, in practice this constraint 
would be further tightened by a more complete calculation. Nonetheless, our analysis 
does not necessarily rule out these models as viable candidates for inflation. Rather, it 
suggests that it will be difficult to construct a working model that results in a detectable 
tensor signal. 
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7. Discussion 

The relativistic DBI brane scenario represents an attractive, string-inspired realisation 
of the inflationary scenario. Recent cosmological data has placed very strong constraints 
on the simplest models based on a single D3-brane. The strength of these constraints 
follows from field-theoretic upper limits on the tensor-scalar ratio, r, which in turn arise 
because the effective DBI action satisfies special algebraic properties. This provides 
motivation for considering generalisations of the scenario, in particular to multi-brane 
configurations. 

In this paper we have identified a phenomenological class of effective actions for 
which the constraints on r are relaxed if significant (and detectable) non-Gaussian 
curvature perturbations are generated during inflation. Included in this class is the 
relativistic limit of the action associated with n coincident branes in the small n limit. 
Moreover we have found that such an effective action for arbitrary, finite n can be 
expressed directly in terms of the corresponding action for the n = 2 model, due to the 
fact that the spin-1/2 representation of SU(2) is actually the fundamental one. This 
allows us to construct models for various values of n using the two-brane action and the 
iteration equations. Physically these brane configurations typically have a smaller sound 
speed than the single brane models due to the different structure of their action. This 
differing structure is also manifest in the non-relativistic limit - since the non-Abelian 
nature of the theory introduces new 'potential terms' that couple to the usual kinetic 
components of the action. In some cases this extra potential could help to further flatten 
the inflaton potential, whereas in other cases it will make it significantly steeper. An 
indepth analysis of slow roll in such models would be welcome. Their backreaction 
is also significantly smaller than other multi-brane configurations, and therefore this 
relaxes the amount of tuning required for the background charge. 

We then proceeded to consider the question of whether the upper limits on r 
could be relaxed to such an extent that a background of primordial gravitational waves 
might be detectable in future CMB experiments. The vast majority of string-inspired 
inflationary models that have been proposed to date generate an unobservable tensor 
background. We found that a detectable signal is possible, in principle, for typical 
string-theoretic parameter values if the number of coincident branes, n, is either 2 or 
3. This is consistent with known F-theory configurations and current WMAP3 limits 
on the non-Gaussianity. Furthermore, we found that the level of non-Gaussianity must 
exceed /nl^ — 70 if such configurations are to generate a detectable tensor signal. This 
is well within the projected sensitivity of the Planck satellite. 

Our analysis invoked an AdS^ x X5 warped throat geometry. However we made 
no assumptions regarding the form of the inflaton potential, other than imposing 
the implicit requirement that the universe underwent a phase of quasi-exponential 
expansion. In this sense, therefore, we have yet to explicitly establish that these 
inflationary models will be able to generate a measurable tensor signal. Nonetheless, 
since such a detection would provide a unique observational window into high 
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energy physics, our results provide strong motivation for considering the cosmological 
consequences of these multi-brane configurations further when specific choices for the 
inflaton potential are made. In particular, it would be interesting to employ the 
techniques developed in [161 EES! [151 El] to identify the ranges of parameter space 
that are consistent with current cosmological observations. It would also be interesting 
to investigate whether the effective action fl29|) with values of m / -1/2 arises in other 
string-inspired settings. 
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Appendix: Exact Solution 



Eq. ( 1251) can be analytically solved in full generality without imposing the limits (12611 
on the derivatives of the kinetic function. This allows us to determine the most general 
class of models where the non-linearity parameter satisfies the condition /nl oc 1/c 2 s at 
leading order. 

In general Eq. (125]) takes the form 

2a 

(2 - a)P, x P, xx + AXP 2 XX = —XP X P >XXX (70) 
and this reduces to 

aQ,x = (6-a)Q 2 + ^^|, (71) 
where Q = P } xx/P,x- Eq. (1711) can be transformed into the linear equation 

'u^*ep>V!!zl (72) 

2a X a 

after the change of variables U = 1/Q and the general solution to Eq. (I72p is given by 

^' xx = \ (7^ 

P x X[/ 2 (0)X(«- 6 )/ 2 «-2] ' { ' 

Integrating a second time implies that 

P x = / 1 (0)(l-/ 2 (0)X- s ) 1/(2s) , (74) 

where s = (a — 6)/ (2a) and we have redefined the arbitrary integration functions fi(4>)- 
Finally Eq. (1741) can be formally integrated in terms of a hypergeometric function 

P = f 1 X 2 F 1 (--,~;l--,hX-') , (75) 
\ s 2s s J 

which represents the most general solution for this class of models. Note that we have 
set the remaining constant of integration to zero to ensure that the kinetic function 
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vanishes in the limit of zero velocity. In fact this expression admits many different 
classes of solution, arising as limits of the expansion of the hypergeometric function. 

The special case of a = 2 (s = —1) implies (after a further redefinition of the 
functions fi{4>)) that 



P = h^\-h^~h (76) 

and this corresponds to the standard DBI action ff23l) [Ml [55]. 

The case a = 18/5 (s = —1/3) can also be expressed in terms of elementary 
functions, again after redefinition of the fi{4>): 



fi 8 - - (J 2 XV3) 

P = — l - 



vnp ~ h - (77) 

Note that this expression appears in a slightly different form to that in (j5"Tl) . However 
in deriving (|5T|) we assumed the relativistic limit, which in turn imposes a non-trivial 
relation between X and (p. Using this, and with a suitable redefinition of the functions, 
we can easily transform the above expression into the required form. 
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